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Abstract 

The presently world largest data sample of 7r° — > 7e+e~ decays containing 
nearly 5xl0 5 events was collected using the WASA detector at COSY. A search 
for a dark photon U produced in the tt° — > — > je + e~ decay from the 
pp ppir reaction was carried out. An upper limit on the square of the 
U — 7 mixing strength parameter e 2 of 5 x 10~ 6 at 90% CL was obtained for 
the mass range 30 MeV < Mjj < 90 McV. This result together with other 
recent experimental limits significantly reduces the M v vs. e 2 parameter space 
preferred by the measured value of the muon anomalous magnetic moment. 
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Figure 1: Fcynman diagrams for a) the lowest order electromagnetic n° — > e+< 
a possible contribution of U vector boson to: b) tt° e+e~7 and c) lepton g 



7 decay and 
2. 



1. Introduction 



Decays of neutral pseudoscalar mesons into a lepton-antilepton pair and 
a photon, P — > are among the tools to search for a new light vector 

boson connected with dark gauge forces [H, 0, H[ . The boson is postulated in 
several theoretical models to be responsible for the annihilation of dark matter 
particles. Such annihilation into an e + e _ pair could explain the positron and/or 
electron excesses observed by PAMELA Q, ATIC @] and H.E.S.S. Q as well 
as the narrow 0.511 MeV 7 ray emission from the galactic bulge observed by 
INTEGRAL Q. 

In one of the simplest scenarios dark matter particles belonging to an addi- 
tional abelian gauge symmetry are added to the Standard Model (SM). The new 
symmetry leaves the SM particles unchanged [1, 0, 0, E3] • The associated gauge 
boson can communicate with the SM through a small mixing in the kinetic term 
of the QED Lagrangian ll| : 



C ■ - -1 pQEDp^ 



(1) 



where e is the mixing parameter. The gauge boson U (or A') is often called 
a dark photon since it can mix with the photon in all processes (examples are 
shown in Figs.[Tp) and c). General theoretical arguments 12, l]| suggest that 
the e parameter must be of the order of 10~ 4 — 10~ 2 and the boson mass Mjj 
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below 1 GeV. This estimate is also supported by the astrophysical observations 
and the constraints imposed by precision measurements such as the anomalous 
magnetic moments (<? — 2) of muon and electron [l4j |. The contribution of the 
U boson to the [g — 2)i (I = e, fi) is given in 14] by: 



A(s-2), = 



2tt 



2m 2 z(l 



2tt 



m 2 (l-z) 2 + M 2 z 



qe 2 2 / m; 
2tt 3 \M L , 



if My < mi 

if Mjj 3> TO; 



(2) 

Investigations of the M\j vs. e 2 parameter space corresponding to the experi- 
mentally preferred (g — 2) M value (shifted by +3 x 10 -9 with respect to the SM 
value 15[| ) are therefore of great importance. 

The U boson is expected to be a narrow resonance. For a U boson with 
mass less than twice the muon mass the total decay width is for all practical 
purposes (neglecting higher-order electric, tiny weak interaction contributions 
from the U boson - Zq coupling, and the decay to light dark scalars and/or 
fermions) given by IS 17 1: 



Trj = Tr 



= iae 2 M v 



1- 



4m 2 



2m z e 



(3) 



where m e is the electron mass. 

Such a light U boson can be directly produced in particle accelerators, see 
e.g. Refs [3 EE El El US IH [H, [3, Q . The idea is to search for narrow 
structures in the invariant mass spectrum of the lepton-antilepton pair. 

The Mjj vs. e 2 region corresponding to the measured (q — 2) M value ±2cr is 
covered by the data from the BABAR [Hj], MAMI Al KLOE-2 [27| and 
APEX [28[ experiments for Mu masses above 100 MeV. On the lower end this 
preferred region is excluded by the (g — 2) e value for Mjj < 30 MeV [29|, |30j. In 
addition, e regions below 10 -6 are excluded by experiments which are sensitive 

,EHH. 
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to lepton pairs from displaced secondary vertices (tjj > 10~ n s) 

Our experiment aims at searching for a short-lived U boson in the ir° Dalitz 
decay, ir° — > e + e _ 7, covering the range preferred by the experimental value of 
(g - 2) M for 30 MeV< M v < 100 MeV. In this region, for e > 10~ 3 the average 
distance passed by a boson emitted from a low energy 7r° decay should be less 
than a millimeter. The best limit from a previous 7r° — > e + e _ 7 experiment with 
the origin of the e + e~ pair close to the production vertex was obtained by the 
SINDRUM collaboration more than twenty years ago (HUff. The SINDRUM 
result is based on a sample of 98400 ir° — > e 
masses above 25 MeV. 



7 decays with e + e invariant 



2. The Experiment 

The WASA detector setup was built and first used at CELSIUS in Uppsala 
and moved to COSY (COolcr SYnchrotron) Julich in the Summer of 2005 [36j]. 
The detector was designed and optimized for studies of rare 7r° meson decays 
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produced in hadronic interactions [37j. It consists of three main components: 
The Forward Detector (FD) - covering scattering angles in the 3° — 18° range 
used for tagging and triggering of meson production, the Central Detector (CD) 
- used for measuring meson decay products, and the pellet target system. The 
target beam consists of 20 - 30 /zm diameter pellets of hydrogen, providing an 
areal target density in the order of 10 15 atoms/cm 2 . The diameter of the pellet 
beam is ~ 3.8 mm. 

The CD surrounds the interaction region and is designed to detect and iden- 
tify photons, electrons, and charged pions. It consists of an inner drift chamber 
(MDC), a superconducting solenoid providing the magnetic field for momentum 
determination, a barrel of thin plastic scintillators (PS) for particle identification 
and triggering, and an electromagnetic calorimeter. The amount of structural 
material is kept to a minimum to reduce the amount of secondary interactions 
outside of the detector sensitive volumes. The beryllium beam pipe (diameter 
6 cm) wall is 1.2 mm thick and the material of the superconducting solenoid 
corresponds to 0.18 radiation lengths. 

The FD allows identification and reconstruction of protons from the pp — > 
ppn° reaction close to threshold. The track coordinates are provided by four 
sets of straw proportional chambers. Kinetic energies are reconstructed using 
the AE information in layers of plastic scintillators of different thickness. In 
addition, the signals are used for triggering. The kinetic energy, T, of the 
protons can be reconstructed with a resolution of u(T) /T ~ 1.5 — 3% for kinetic 
energies below 400 MeV. 

The results presented here are based on data collected during one-week 
WASA-at-COSY run carried out in 2010. The tt° mesons were produced in 
proton-proton interactions at a kinetic beam energy of 550 MeV. The beam en- 
ergy corresponds to the center-of-mass excess energy of 122 MeV with respect to 
ppir° threshold (i.e. below two pion production thresholds) with a cross section 
of 1.12 mb [381 ] . The maximum scattering angle of the outgoing protons for the 
reaction is 45° . For detection and for triggering purposes the phase space of the 
pp — > ppn reaction can be divided into three regions: 

1. Both protons are measured in the FD. This corresponds to a geometrical 
acceptance of 19%. 

2. One proton is measured in the FD and one in the forward part of the PS 
(scattering angles 20°— 40°). This corresponds to a geometrical acceptance 
of 42%. 

3. Both protons are registered in the PS. This corresponds to a geometrical 
acceptance of 21%. 

Case (1) allows the definition of the most selective trigger condition and the 
best resolution in the missing mass with respect to the two protons. Therefore, 
the main trigger for the experiment required two tracks in the FD. The protons 
from the pp — > ppn reaction have a maximum kinetic energy of 350 MeV and 
are mostly stopped in the FD. This allows the inclusion of a veto from a thin 
plastic detector layer placed at the far end of the FD into the trigger condition. 
In addition, two hits in the central part of the PS (scattering angles 45° — 135°) 
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were required, aiming to select the electron-positron pair. An additional, scaled 
down, trigger based on case (2) was used in parallel. The WASA-at-COSY data 
acquisition system allowed the collection of more than 10 4 events per second 
and the luminosity was set to optimize the conditions for the main trigger. The 
integrated luminosity of the run was about 0.55 pb _1 . 

The data quality is illustrated by analysis of the main trigger data sample 
and requesting in the analysis two identified (using AE/ AE method) FD proton 
tracks. An electron positron pair is selected by requiring two oppositely curved 
tracks in the MDC with scattering angles between 40° and 140°. A photon 
hit cluster in the calorimeter with an energy deposit above 20 MeV is also 
requested. The missing mass squared with respect to two protons (MM 2 (pp)) 
for the above selection is shown in Fig. [2ji. In addition to the pp —> ppir reaction 
signal one sees also a contribution due to random coincidences of pp — > pp and 
pp — > pnir + reactions. This background is effectively suppressed by including 
electron and positron identification using the reconstructed momentum and the 
energy deposit in the calorimeter. The corresponding MM(pp) plot after this 
cut is shown in Fig.[5j3. The 7r° — > e + e~7 decay is independently identified from 
the invariant mass of the decay products IM(e + e~j) (calculated assuming the 
tracks originate at the beam target crossing) shown in Fig. [5J:. The data are 
well described by a simulation of pp —> ppn with ir° — > e + e~7 and 7r° — > 77 
decays, where in the latter case one of the two photons converts in the beryllium 
beam tube. 

For Monte Carlo simulations, angular distributions for the pp — > ppir reac- 



tion from 38[ were used in the event generation. The tt -> e + e 7 decay is 



generated using the lowest order QED matrix element squared: 

\A\ 2 = r 77 16^M^l (l - (l + cos 2 r + lpisin 2 r^) \F(q 2 )\ 2 (4) 

where 9* is the angle of e + in the dilepton rest frame with respect to the dilepton 
momentum in the overall ir° decay system, M and m e are 7r° and e ± masses 
respectively, T 77 is the partial 7r° — > 77 decay width, and F(q 2 ) (with q 2 the 
squared momentum transfer of the off-shell photon) is the ir° transition form 
factor. The form factor close to q 2 — is parametrized as: F(q 2 ) = 1 + 
aq 2 /M 2 . The value of the dimensionless linear coefficient a is 0.032 ±0.004 (39[. 
Consequently, for the range of q considered in this work the form factor can be 
safely approximated by one. 

The matrix element from Eqn. Q leads to the following unperturbed dT/dq 
distribution [4(J for the standard lowest order electromagnetic decay tt° — > 
e + e _ 7 of Fig. Q}i: 
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c ) IM(eVy) [GeV] 



Figure 2: Control data sample with two protons, an e+e~ pair and a photon, a) Distribution of 
the missing mass squared with respect to the two protons registered in the FD before electron 
identification. Experimental data (black points); simulations: n° — > e+e~ 7 and w° — ► 77 
(broken line), random coincidences of two events (dotted line), and the sum (solid line), b) 
Distribution of MM(j>p) after electron identification: experimental data (black points) and 
sum of Monte Carlo simulations (solid line), c) The reconstructed invariant mass of the e+e~7 
system after particle identification cut. 
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b ) IM(e + e ) [GeV] 

Figure 3: The reconstructed e+e — invariant mass q = 7M(e+e — ): a) before and b) after 
the cuts for reducing the conversion background. The experimental data are denoted by 
black points. Results of simulations for n° — > 77 (blue line) and iz° — > e+e — 7 (green line) 
decays are normalized according to the known branching ratios. The normalization of random 
coincidences (dotted line) was fitted in order to reproduce the IM(e + e~) > 150 MeV range. 
The sum of all simulated contributions is given by the red line. 



3. Data analysis 

The first stage of data analysis is to extract a clean signal of ir° — > e + e~j 
decays. The results shown in the previous section suggest that at this en- 
ergy electron-positron pairs come nearly exclusively from the 7r° meson decays. 
Therefore, in order to maximize the yield of the tt° — > je + e~ events we use 
an inclusive data sample requesting events with (i) at least one proton iden- 
tified in the FD, (ii) an e + e~ pair identified in the CD. There is no request 
of an additional photon cluster and we have included events from both trig- 
gers corresponding to phase space regions (1) and (2). The distribution of the 
reconstructed invariant mass of the electron-positron pair, q = IM(e + e~), is 
shown in Fig. [3^. This spectrum is well described by the sum of 7r° — > e + e~7 
and 7T° — > 77 (with photon conversion). The data sample contains 1.8 xlO 6 
reconstructed events. 

The 7T° — > 77 events are efficiently removed by a condition on the recon- 
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10 15 20 25 30 35 40 45 50 

R [mm] 



Figure 4: Distribution of the distance R between the COSY beam axis and the reconstructed 
point of closest approach of e + and e~ tracks: experimental data (black crosses); simula- 
tions for 7T° — ¥ 77 (blue line), the ir° — ► e+e — 7 decay (green line), and the sum of the two 
contributions (red line). 



structed position of the e + e~ vertex. Fig. 0] shows the distance (R) of the re- 
constructed vertex from the COSY beam axis. The contributions of the ir° — > 77 
and 7T° — »■ e + e~7 decays, simulated according to the known branching ratios, 
are in very good agreement with the observed distribution and they are well 
separated. In order to further reduce the external conversion background one 
uses the invariant mass of the e + e~ calculated from the momentum directions 
at the points where the tracks intersect the beam tube, 1Mb, shown in Fig. [SJ 
The selection cut is performed in the 1Mb vs. R plane (Fig. [5]). The cut removes 
98% of the 7T° — > 77 events which contribute to IM(e + e~) distribution due to 
conversion. 

The finally reconstructed dN/dq distribution, containing nearly 5xl0 5 en- 
tries, is shown in Fig. [5Jd. It is well described by the simulations of the 7r° — > 
e + e~7 decay channel alone with a very small (approx. 3000 events) admixture 
of background from the ir° — > 77 decay. The data in this work represent the 
world largest data sample of ir° — > e + e~j events, which is almost an order of 
magnitude larger than the sample used for the previously published results from 



mip le 1 
the SINDRUM experiment 



3.1. Upper limit for the BR(n° -!• -f(U -> e+e )) 

A distinctive feature of the expected signal of the decay 7r° — >• j(U — > e + e~) 
(Fig. lb) is the appearance of a narrow peak (the width being given by the 
detector resolution) in the invariant mass distribution of the electron positron 
pair at the U boson mass. The electrodynamics process 7r° — > 7*7 — > e + e~j 
(Fig. la) both represents the irreducible background and is used for normaliza- 
tion. Due to the expected small decay width of the U boson the interference 
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R [mm] 



Figure 5: Correlation between R and IMf, variables for the experimental data. The selection 
cut is shown by the diagonal line. The events below the line mainly come from photon 
conversions in the beam pipe. 
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Figure 6: Difference between the reconstructed e + e — invariant mass distribution and the sum 
of all simulated contributions (black points). The resolution and sensitivity are illustrated by 
the superimposed red histograms. They represent the signals expected for the ir° — > U-f — ¥ 
e + e — 7 process with U boson masses of Mu = 30, 50, 70 and 90 MeV and BR(n° — » U^/) = 
10 -4 (the corresponding e values are: 0.0077, 0.0088, 0.0113, and 0.0169 respectively). 
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term is negligible and the signal from the U boson can be tested by constructing 
an incoherent sum of the two contributions. 

The experimental data are described well by the simulation based on Eqn. Q 
alone as shown in Fig. [3] The difference between reconstructed experimental 
q distribution and the sum of all simulated contributions is given in Fig. [51 
In addition there are superimposed five example distributions corresponding to 
the 7T° — > — > e + e~7 process for U boson masses of 30, 50, 70 and 90 MeV 
respectively, assuming BR(ir° — > Uj) = 10~ 4 . The plots illustrate both the 
resolution and the efficiency expected for the signal. 

For a given value of U boson mass corresponding to the range of k th bin of the 
invariant mass spectrum (qk < Mjj < qk+Aq, with Aq = 1 MeV the width of the 
histogram bin) the number of events in the i th bin of the reconstructed electron- 
positron invariant mass distribution, iVj, can be described in the following form: 

1 



Ni/N Tot = - V" SijTjjVj + Sikiikfi (6) 



The first term in the Eqn. ([5]) represents the contribution from the Dalitz decay 
and the second term from the hypothetical 7r° — > j(U — > e + e~) decay chain. 
Indices j and k label the true, unperturbed distributions and i the reconstructed 
q histogram. Nrot is total number of produced 7r° mesons, 1/r is the 7r° life 
time and rjj is the efficiency. Sy is the normalized smearing matrix (for each 
j: '^ li Sij — 1), Vj is the unperturbed dT/dq distribution for the 7r° — ^ e + e~7 
decay (Eqn. ([5])) integrated over bin j: 

+Aq dT 

=-dq, (7) 

and (3 is BR(tt° — > j(U —> e + e~)). The efficiency and the smearing matrix 
was obtained from the detector simulation. The U boson decay mechanism in 
diagram Fig.[T]3 implies that the efficiencies as a function of cos 6* are identical 
to the ones of the ir° — > e + e~7 decay with q = Mjj. Note that for the quoted 
values of the branching ratios the intrinsic width (3) of the U boson would be 
in the eV range and thus very much smaller than the experimental bin size. 

The upper limits for the U boson branching ratios, /?, as a function of Mjj 
were obtained by repeating for all bins (index fc in Eqn. ([5])), corresponding to 
the 30 MeV < M v < 100 MeV range, the least square fits of Eqn. © to the 
experimental q distribution. Fig.[7]shows the 90% CL upper limits obtained for 
the branching ratio of ir° — > -f(U — > e + e~) decay as a function of the assumed 
value of Mjj- This result is compared to that obtained from the SINDRUM 
data [jyj. 



The branching ratio of tt° — > jU is related to e by [42, ll 



} =2e 2 \F(M*)\ 2 l--r^ . (8) 



r(7r ^ 77 ) y M 2 

The resulting upper limits for the e 2 parameter is shown in Fig.[8]and compared 
with other experiments. 
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Figure 7: A 90% C.L. upper limit for the BR(jt° — > 7I/) from this paper (solid line) compared 
to the result of the SINDRUM experiment [3 (dotted line) . 

The recent limits for the electron g — 2 are taken from Refs [3(], [29| . Our 
upper limit is compatible with the combined KLOE limits [43j using a disparate 
experimental setup and different meson decay as source of e + e~ pairs. Together 
they significantly reduce the parameter space for mass and mixing strength of a 
hypothetical dark photon U, if the latter is assumed to account for the present 
experimental value of the muon anomalous magnetic moment. The experiment 
presented in the paper if repeated with an order of magnitude larger statistics 
would cover the remaining part of this region of interest. The collected data 
can also be used for to determine the ir° transition form factor. 
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Figure 8: Summary of the 90% CL upper limits for the mixing parameter e 2 from WASA- 
at-COSY (red solid line) compared to SINDRUM 7r° -> e+e~7 [34| (dotted line) and recent 
combined KLOE <j> —>■ r)e~^e~ [43| (dashed dotted) upper limits. The long respectively short 
dashed lines (and the corresponding hatched areas) are the upper limits derived from the 
muon and the electron g — 2 [29j . In addition the gray area represents the ±2<r preferred band 
around the present value of the muon g — 2. 
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